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1. INTRODUCTION 

The Wiener model is a series connection of a linear dynamic bloc and a memoryless nonlinearity 
(Figure 1). Nonlinear systems structured in block are widely studied in the last decades [1-10]. When both 
parts are parametric, the identification problem has been dealt with using stochastic methods, e.g. [1, 2] as 
well as deterministic methods, e.g. [11, 12]. The stochastic methods enjoy local or global convergence 
properties under various assumptions e.g. the system inputs should be persistently exciting (PE) or Gaussian 
and the system nonlinearity is invertible. The last limitation has recently been overcome in [12]. Multi-stage 
methods, involving two or several stages, have been proposed in e.g. [13, 14] and their consistency was 
ensured if the inputs are Gaussian and the nonlinearity 1s odd. Deterministic parameter identification methods 
consist in reformulating the problem as an optimization task, generally coped with using various relaxation 
techniques. Then, local convergence properties ensured in presence of PE inputs. In Nonparametric Wiener 
systems (where none of the linear subsystem or the nonlinear element assumes a priori known structure) have 
been dealt with using both stochastic and frequency methods. In the stochastic methods e.g. [15], 
the nonlinearity is generally determined using variants of the kernel regression estimation technique while the 
(unknown) coefficients of a FIR/IR approximation of the linear part are estimated using cross-correlation 
analysis. Several assumptions are needed e.g. Gaussian inputs, FIR linear dynamics, Lipschitzian 
nonlinearity. In frequency methods, the linear subsystem frequency response and the nonlinearity map are 
determined in two or several stages, e.g. [16, 17]. The frequency identification of nonlinear system has been 
widely dealt by Brouri for several structures and models [16-19]. Presently, the considered nonlinear system 
is characterized by nonparametric linear block and can be of unknown structure. Furthermore, the system 
nonlinearity is of hard type (Figure 2). This latter may have several phenomena; e.g. saturation, preload, dead 
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zone. Presently, the system nonlinearity considered in Wiener models is very difficult to be modeled or 
approximated by polynomial function or any other basis functions. 

The aim in this work is to develop an identification method allowing to estimate the hard system 
nonlinearity, as well as an accurate estimate of linear block parameters. The proposed identification approach 
can be developed in two stages. In the first stage, using simple constant inputs an accurate estimate of 
the nonlinearity functions f(.) can be given. In the second stage, the frequency modulus gain |G(jæ)| and 
the argument (the phase) 4G (jw) of G (jæ), for a set of frequencies w E {Wy , Wo, ..., Wm}. 

The outline of the remaining part of this paper consists of 3 sections. The system nonlinearity 
identification is coped with in Section 2. The determination of the linear subsystem frequency response is 
investigated in Section 3. Section 4 is devoted to the conclusion. 





Figure 1. Wiener Model structure Figure 2. Nonlinearity hard with preload and dead zone 


2. IDENTIFICATION OF SYSTEM NONLINEARITY 

Standard Wiener models consist of a linear element of transfer function G(s) followed in series by 
a static nonlinear function f(.). The Wiener systems considered in this work is characterized by set of 
internal signals v(t) and w(t), and €(t) (Figure 1), which are inaccessible to measurement. The external 
signals u(t) and y(t) respectively denote the control input and the measured output. The inner signal (t) 
accounts for measurement noise and other modeling effects. It is supposed to be zero-mean ergodic and 
uncorrelated with the control input u(t). Accordingly, let g(t) denotes the impulse response of the linear 
block of transfer function G(s). 

More specifically, the Wiener system under study is analytically described by the following 
equations: 


v(t) = g(t) * u(t) (1) 
w(t) = f(v(t)) (2) 
y(t) = w(t) + E(t) (3) 


where g(t) = L(G (s)) and L~* designates the inverse Laplace transform. The symbol * in (1) refers to 
the convolution product. 

Presently, the considered linear block can be of unknown structure, but it is just supposed to be 
asymptotically stable. This assumption makes possible open-loop system identification. In this section, an 
identification schema is proposed to estimate the hard system nonlinearity (Figure 2). Then, the Wiener 
system is excited by a set of constant inputs: 


u(t) €{U;;j =1..n} (4) 


where the number of input excitations n in (1) is arbitrarily chosen by the user. 

In this problem identification, the complexity also lies in the fact that the internal signal v(t) is not 
uniquely defined from an input-output viewpoint (see e.g. [16-19]). Then, without loss of generality, one can 
assume that the linear block to be identified is characterized by unit static gain: 


G(0) = 1 (5) 


Indeed, if the couple (G(s), f (.)) is solution of this nonlinear identification problem (1)-(3), then any 
nonlinear system having the following linear subsystem block: 
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co% (6a) 


and the following system linearity: 
f°) = fav) 


is also solution of the above nonlinear identification problem (1)-(3), where is any nonzero real number. 
Accordingly, the judicious choice of the distribution factor in the present paper is: 


A = G(0) (6c) 


Finally, one immediately gets from (6a) and (6c) that, the linear block G*(s) to be identified is of unit 
Static gain: 


G*(0) =1 (6d) 
On the other hand, it is readily seen using (1) and (4)-(5) that: 
v(t) € {U;;j =1..n} (7) 


Equation (7) means that, for any constant input U; the steady state of inner signal v(t) is also constant and 
converges towards U; (i.e. of known values). Accordingly, one can immediately gets from (2) and (7) that 
the undisturbed system output w(t) converges towards set constant values. Specifically, one has: 


w(t) €{f(U;);j =1..n} (8) 


Finally, it follows from (3) and (8) that, for any input U;, the system output y(t) converges to a constant up 
to noise: 


y(t) EFU) + FO; =1...n} (9) 


Then, using the fact that, the extra-input (noise) signal €(t) is zero-mean, an accurate estimate of f (U;) 
(for j = 1... n) can be obtained by averaging the system output y(t). To this end, let Wi = 1... n) denotes 
the estimate of w(t) for u(t) = U;. Then, the estimator of w(t) for any j € {1, ..., n} is as follows: 


N 
so a 
W; = W220 (10) 


Accordingly, it readily follows from (3), (9) and (10) that: 
j N 
W, = f(Uj))+—) EC) (11) 
t=1 
One can see that the last term in (11) vanishes with probability equal to 1. Accordingly, the estimate W; can 
be expressed as (w.p.1): that: 
W, ~ f(U;) forj=1..n (12) 
For more information, see e.g. [16-19]. Accordingly, a set of points belonging to nonlinearity function f (.) 


can be determined. 


3. RESEARCH METHOD 

The aim in this section is to develop an identification solution allowing estimates of linear block 
parameters. Presently, the proposed Wiener model is characterized by nonparametric linear element. This 
latter can be of unknown structure. Accordingly, the objective of this identification method provides 
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the frequency complex gain of linear subsystem, i.e. for any frequency w E { w4, ... Wm} we seek an estimate 
of the modulus gain |G(jw)| and the phase p(w) = 4G (jw) = arg(G (jw)). 
To do this, the Wiener model (Figure 1) is excited by a sine input signal: 


u(t) = ug + Ucos(wt), with w E { %4, ... Wm} (13) 


Then, it follows from (1), (5) and (13) that, the steady state of inner signal v(t) is also sine signal and can be 
expressed as follows: 


v(t) = uo + U|IGGw)| cos(wt + p(w)) (14) 


On the other hand, recalling that the system nonlinearity f(.) is of hard function and is already determined 
(Section 2). Then, the choice of ugo and U can be done such that the signal v(t) spans only one segment of 
f(.). In that case, the steady-state of undisturbed system output w(t) is also sine signal. Specifically, w(t) 
can be expressed as: 


w(t) = uo + b+ au|GGw)| cos(wt + p(w)) (15) 


where a and b are the segment parameters. 

A judicious choice of ugo and U can be made practically by observing the system output y(t). Then, 
adjusting the offset and the amplitude until y(t) becomes close to a sinusoidal signal up to noise. 
Accordingly, it readily follows from the (3) and (15) that, the system output y(t) can be expressed as: 


y(t) = uo + b + aU|GGw)| cos(wt + p(w)) + E(t) (16) 


Furthermore, the question in this stage is how to estimate the inner signal w(t) using input and 
output system? The answer to this question can be found in [20]. Know that the input signal u(t) and 
the system output y(t) are T — periodic (of same period T). Then, an accurate estimate W(t) of w(t) can be 
obtained by trans-averaging y(t) (see [20]): 


M 
1 
w(t) = =) y(t +17) (17) 
l=1 


Furthermore, it follows from (16) and (17) that, the output estimate W(t) can be expressed as: 


M 
w(t) = uo + b+ aU|GGw)| cos(wt + p(w)) + =). E(t + IT) (18) 


l=1 


Bear in mind that, the signal €(t) is a white noise process with zero mean. Then, the expression of W(t) in 
(18) can be reduced to the following equation as: 


w(t) = uo + b+ aU|GGw)| cos(wt -+ p(w)) (19) 


Equation (19) shows that the estimate W(t) of internal signal w(t) converges to its true value (w.p.1). 
Specifically, one has: 


w(t) =w(t) w.p.l (20) 


Let S(w) and @(w) respectively denote the amplitude and phase of the sinusoidal term of W(t). 
Then, the frequency gain estimate (modulus gain |G(jw)| and phase @(w) = ZG(jw)) can be given as: 


|GGw)| = 2 (21a) 
@(w) = &(w) (modulo 27) (21b) 


This experience can be repeated for any other frequency w E { w1, ... Wy}. 
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4. SIMULATION 
The Wiener systems considered in simulation is described linear subsystem of transfer function: 


EPE 0.1 32 
ETE) a 


Then, the system nonlinearity f (.) considered in simulation shown by Figure 3. Accordingly, to overcome 
the multiplicity of solutions problem seen in section 2, we have shown that it is judicious to focus to 
nonlinear system described by (6a-c). Then, the nonlinear system to be identified is analytically described by 
linear element of following transfer function: 


aie 0.04 D 
‘= Grane +04) 


Furthermore, it is readily seen from (23) that the linear element of Wiener nonlinear system to be 
identified satisfies the property (5), i.e. G*(s) is of unit static gain (G*(0) = 1). In the first stage, the Wiener 
system is excited by a set of constant inputs (Figure 4). The resulting system output y(t) and the undisturbed 
output estimates W, are shown by Figure 5. For convenience, the system nonlinearity f(.) and the set of 
estimated points (U;, W;) (Vj =1...9) are plotted in Figure 6. Finally, the obtained simulation values and 
given in Table 1 shown that the estimated points are very close to their true values. 


The hard nonlinearity 
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Figure 3. The nonlinearity considered in simulation Figure 4. The input signal u(t) 
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Figure 5. The output y(t) and the estimates W; Figure 6. The NL f (.) and the estimates (U;, W;) 


Table 1. Nonlinear points estimates f (U;) 


J 4 5 7 
U; -0.25 0 0,5 
f(U;) -0.87 0 0,25 
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Î(U;) -0.88 0.02 0,27 

In the second stage (see section 3), the Wiener nonlinear system is excited by the sine signal (13). 
For an arbitrarily offset uo and amplitude U, the system output y(t) over one period is shown by Figure 7. 
Then, it is clear that the signal y(t) is different from a sinusoidal signal. 

Accordingly, by adjusting the values of wg and U in (13) and using the nonlinearity estimate. 
Furthermore, using the system nonlinearity estimate (Figure 6), useful values for uo and U in (13) can be 
easily done. Then, another example for wy = 0.25 and U = 0.5 is displayed by Figure 8. 

This result showed that the system output is not yet a sine signal. Indeed, the filtered version using 
the estimator (17) is presented in Figure 9. It is seen that the signal y(t) is different from a sine signal. Then, 
adjusting the parameters of input signal u(t) (i.e. the values of offset uọ and amplitude U) until the system 
output y(t) becomes sine signal (up to noise). The observed system output y(t) in this case is illustrated by 
Figure 10. 


The system output y(t) The system output y(t) 
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Figure 7. The system output y(t) over one period Figure 8. The output y(t) over one period 


The filtered output version 


The system output signal y(t) 
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Figure 9. The filtered output W(t) over one period Figure 10. The output y(t) over one period 


This result shows that the signal y(t) is close to a sine signal (up to noise). The collected sample of 
y(t) is used to generate the undisturbed output estimate W(t) using the estimator (17). This latter is plotted 
by Figure 11. It is remarkable that W(t) is close to a sine signal. Then, for the frequency w = 0.01(rd/s), 
the obtained amplitude and phase of W(t) are as follows: 
(w) = 1.01 a 


ĉâ(w) = 3.04 (rd) (24b) 
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Then, based upon the system nonlinearity identification, the frequency complex gain (modulus gain |GGo)| 
and phase @(w) = ZG(jw)) of linear block can be given using (21a-b). 

Bear in mind that â ~ —0.99, then it readily follows from (21a-b) and (24a-b) that, the parameters 
estimate (modulus gain and phase) of frequency complex gain for w = 0.01(rd/s) are (respectively): 


|GGjw)| = = = 1.01 (25a) 
@(w) = 3.04 — r ~ —0.1 (rd) (25b) 


Repeat this method for other frequencies w E { w1, ... Wy }.S 

It follows (13a-b) and using the amplitude and phase of W(t), can be obtained for any frequency w. 
The obtained estimate values of the modulus gain |G(jw)| and phase @(w) or three frequencies w € 
{ W1,... Wy} are given in Table 2. Finally, the obtained results given in Table 2 show that, the complex gain 
parameters estimate are very close to their true values. 


The filtered output W(t) 


Tt t T Tt 


Table 2. Frequency gain estimates G(jw;) 





k l 2 3 
0.5 W, (rd/s) 0.01 0.05 0.1 
~ (Wx) (rd) -0.12 -0.59 -1.03 
0 @ (Wx) (rd) -0.1 -0.62 -1.06 
IGGax)| 0.99 0.89 0.68 
05 |GUox)| 1.01 0.87 0.71 
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Figure 11. The undisturbed output estimate W(t) 
over one period 


5. CONCLUSION 

In this paper, a new frequency identification method is proposed to deal with Wiener model. 
Presently, the identification problem is addressed in presence of hard nonlinearity and linear dynamic 
element of unknown structure. In the first stage, an identification method is presented to estimate the hard 
nonlinearity. Then, a frequency-domain identification solution is developed allowing the determination of 
frequency complex gain (i.e. modulus gains estimate, as well as the phase) estimate of linear block.Presently, 
the linear element is nonparametric and can be of unknown structure. Furthermore, the system nonlinearity 
can be discontinuous. 
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